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Abstract. A new method is described that permits quickly and easily, a 2-dimensional search for TeV 7-ray 
sources over large fields of view (~ 6°) with instruments utilising the imaging atmospheric Cerenkov technique. 
It employs as a background estimate, events normally rejected according to a cosmic-ray background rejection 
criterion based on image shape, but with reconstructed directions overlapping the source of interest. This so-called 
template background model is demonstrated using example data taken with the stereoscopic HEGRA System of 
Cerenkov Telescopes. Discussion includes comparisons with a conventional background estimate and limitations 
of the model. The template model is well suited to the search for point-like, moderately extended sources and 
combinations thereof, and compensates well for localised systematic changes in cosmic-ray background response. 
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1. Motivation & Introduction 

The search for new astrophysical sources of TeV 7-ray 
emission today is well motivated by the multi-wavelength 
picture of candidate sources, and also by the high perfor- 
mance offered by present and future instruments. Many 
of the >150 unidentified EGRET sources (MeV - GeV 
energies) have positional uncertainties approaching sev- 
eral degrees (Hartman et al. H999|) . presenting exciting op- 
portunities for TeV instruments to identify possible coun- 
terparts over the coming years. The ground-based imag- 
ing atmospheric Cerenkov technique presently offers on- 
axis flux sensitivities better than (for energies > 1 TeV) 
~ 10~ 12 erg cm~ 2 s _1 for 50 hrs exposure, and can achieve 
arc-minute scale angular resolution over fields of view ap- 
proaching 5° diameter (see e.g. Konopelko et al. 2002). 
Rather wide surveys for new TeV sources are therefore 
possible even from singly pointed observations. Techniques 
to generate 2D skymaps of event excess significance for 
such surveys have therefore become more important in 
recent years. Critically important, for such skymap gen- 
eration is an estimate of the cosmic ray background over 
the field of view (FoV). Many expected, and known sites 
of TeV 7-ray production, such as supernova remnants, 
pulsars/plcrions, microquasars, and nearby extragalactic 
sources will likely present complicated morphology for fu- 
ture instruments as sensitivities improve. Particularly in 
cases where combinations of such sources might be ex- 
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pected in the FoV, it is vital that the background estimate 
is not contaminated by other TeV sources and various sys- 
tematic biases are well understood. 

Described here is the template background model (an 
earlier version is described in Rowell 2000), designed to 
provide a background estimate for sources of interest at 
all positions in the FoV. The model is demonstrated on 
archival data from the HEGRA System of Cerenkov tele- 
scopes (HEGRA CT-System). For detailed descriptions of 
this instrument and its performance see Daum et al.[1997 
and Puhlhofer et al. 2003 Comparisons are made between 
results of the template model and those taken from a con- 
ventional background model presently in use in HEGRA 
CT-System data analysis. Some features and limitations 
of the template model are described. 



2. Conventional Background Estimates and 2D 
Skymap Generation 

Present ground-based instruments utilising the imaging 
atmospheric Cerenkov technique (100 GeV to ~30 TeV), 
must detect 7-ray initiated events against a background 
of vastly outnumbering isotropic cosmic-ray (CR) initi- 
ated events. First, a set of so-called 7-ray cuts are em- 
ployed to preferentially select events conforming to a 7- 
ray hypothesis. These are a-priori chosen, and based pri- 
marily on Cerenkov image shape parameters, for example 
width, length (Hillas I1985|l . and mean- scaled-width (w ) 
fKonopelko 119951 Daum et al. 1997). Second, event ar- 
rival directions are also reconstructed based on the ori- 
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Fig. 1. Illustration of some conventional background es- 
timation methods. (1) Separate OFF source observations 
are taken at matching zenith angles to those of ON source 
data. (2) Discrete background regions within the same 
FoV at a common source-to-tracking distance so as to 
match background response. (3) A continuous version of 
(2), ring or ring segments may be used instead. Regions 
adjacent to the source (shaded) may be preferred to avoid 
global systematic changes in response. (4) A ring region 
completely surrounding the source. Methods 2,3 & 4 are 
termed displacement background models. Method 4 is 
used when the source is too close to the tracking position 
to implement methods 2 & 3, or over the entire FoV. 



entation of recorded Cerenkov images (see eg. Buckley et 
al.^998, Hofmann et al. 1999) and a so-called directional 
cut made on the distance 6 (or analogous parameters) be- 
tween the reconstructed and assumed arrival directions. 
To assess the statistical significance of any region in the 
FoV, an estimate of the CR background surviving both 7- 
ray and directional cuts must be made. This background 
estimate can be derived from separate observations of dif- 
ferent OFF tracking positions (ON/OFF tracking), or con- 
trol regions displaced within the same FoV containing the 
source region. The latter types of background models are 
termed here displacement backgrounds, and a number of 
different displacement background geometries are avail- 
able. Background regions are chosen so that their CR re- 
sponse matches as closely as possible, that of the source 
region (see Fig. In ON/OFF tracking this is achieved 
by choosing background regions with the same relative 
displacement from respective tracking positions as that 
for the source. For displacement backgrounds, such re- 
gions are mirrored through the tracking position (meth- 
ods 2 & 3), except for sources very close to the track- 



ing position in which a region surrounding completely the 
source can be implemented (method 4). In this latter case 
an extra correction is applied to account for differences 
in CR response. Provided the FoV is large compared to 
the angular resolution and source size, and global system- 
atic changes in CR response are minimal or can be al- 
leviated, displacement backgrounds may be implemented. 
Displacement backgrounds are favoured over the ON/OFF 
method since they do not require extra time taking ded- 
icated OFF source observations. A recent modification of 
the ON/OFF method, however, is proposed by Petry|2003 
The large cameras (~4.3° diameter) and stereoscopic trig- 
ger employed by the HEGRA CT-System achieves a quite 
flat FoV (FWHM ~ 3°) that is amenable to displace- 
ment background models over the entire FoV. In this pa- 
per, using example data from the HEGRA CT-System, 
results from a combination of displacement methods 3 
and 4 (termed here the displacement background model) 
are used to compare with the template model. For source 
to tracking distances larger than ~ 0.4°, a ring segment 
or arc (method 3 ) adjacent to the source region (within 
position angles 135° to 165°) is used, otherwise, method 
4 is chosen. Ring background regions rather adjacent to 
the source are often preferred to avoid global systematic 
changes (of ^few percent) in CR response across the FoV. 

For a given region in the FoV, the source events s are 
derived from within a nominal source region according to 
some directional cut 6 < (9 cu t and 7-ray cuts on image 
shape parameters. The number of background events b 
are derived according to the above-described displacement 
background methods, after also applying the same 7-ray 
cuts. The statistical significance S of excess counts s — ab 
is then estimated from Eq. (17) of Li & Ma (|1983f> : 



S = V2 s In 



s(l + a) 



a(s + b) 



bin 



6(1 + a) 



(1) 



where a normalisation a accounts for the solid angle 
ratio between the source and background regions when 
displaced backgrounds are used, or different observation 
times if ON/OFF tracking is used. Under the displaced 
background models applied to HEGRA CT-System data, 
typical normalisation values of 0.05 to 0.3 are achieved, 
thereby considerably improving the background estimate 
compared to an a = 1.0 situation usually achieved in 
ON/OFF tracking. 

In 2D skymap generation, the above procedure is 
repeated over the FoV at a series of grid positions. 
Directional parameters such as 9 and those necessary to 
form the background estimate are then re-calculated with 
respect to each grid position, and excesses and signifi- 
cances are calculated accordingly. 

3. The Template Background Model 

The template background model invokes a different philos- 
ophy to that of displacement background models. Instead 
of selecting spatially different regions in the FoV for back- 
ground estimates, template model background events are 
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Fig. 2. Normalised distributions of w for 7-ray and 
cosmic-ray (CR) events, based on Monte Carlo simulations 
(solid histograms) and real data (error bars). Adapted 
from Konopelko 2002 



comprised of those with their arrival directions recon- 
structed in the same region as that for the source, yet 
separated in image shape parameter space. Such events 
can form a suitable template of response for 7-ray-like CR 
events over the FoV. 

As applied to data of the HEGRA CT-System, the 
shape parameter mean- scaled-width (w ) is used to con- 
struct the template background estimate, w is the image 
scaled-width, averaged over all images accepted for anal- 
ysis. The scaled-width is the image width scaled accord- 
ing to a 7-ray hypothesis. This scaling is dependent upon 
the total photoelectron yield of the image, the zenith an- 
gle of observations, and the impact parameter or distance 
to the air-shower core. By design (see Fig. FJ| , true 7-ray 
events have w values centred at unity, leaving the generally 
broader CR events to fill higher values. In the search for 
weak, background-dominated sources in CT-System data, 
the cut w <1.1 has been found optimal. Template back- 
ground events should be derived from a w regime con- 
taining little or no true 7-rays, i.e. w >1.1. Experience 
has found that events within the regime 1.3 <w < 1.5 
work well, although regimes containing higher w values 
have proved adequate. To apply the template model prac- 
tically however, two dominant systematic effects must be 
corrected: 

— Radial Correction: Differences in the radial CR ac- 
ceptance between the two w regimes. 

— Zenith-Correlated Correction: Differences in CR 
acceptance between the two w regimes which correlate 
with the zenith angle of reconstructed event directions. 

For a nominal region in the FoV defined by a cut 9 < 6 cut , 
7-ray-like events s are selected according to w <1.1. 
Similarly after applying the same cut on 9, template back- 



ground events b are defined as those satisfying the con- 
dition 1.3 <w < 1.5. In summary template background 
events are assigned a weight uio 



1.0 



(w < 1.1) 



(2) 
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rack,: 
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track, i/ 



/(Azi) (1.3 < w < 1.5) (3) 



for events i in the source region. The correction terms 
applied to b can be described with suitable parametrisa- 
tions as follows. The radial acceptance of the 7-ray-like 
(w <1.1) and template (1.3 <w < 1.5) regimes over the 
entire FoV are sufficiently fitted out to radius R, with 8th- 
order polynomials, p s and pb respectively and normalised. 
Both polynomials are functions of the event distance to 
tracking position #track- A strong increase or pile- up of 
reconstructed event directions due to camera edge effects 
is normally seen at the FoV edge. Ignoring Cerenkov im- 
ages whose centroid distance from the camera centre dis 
exceeds 1.7° effectively suppresses this pile-up, and allows 
fitting of the FoV response to the physical camera lim- 
its. The choice of R is limited by event statistics, and 
practically, can reach values up to 3°. Fitting such poly- 
nomials assumes that the radial response is azimuthally 
symmetric. However, a second systematic effect, a gradi- 
ent /(Az) correlating with the difference in zenith angle 
Az between the event and tracking position, is also gen- 
erally present. This can be described by a linear function 
vs. Az, of the ratio of FoV responses between the 7-ray- 
like and template w regimes. This linear fit is made after 
first correcting for the radial profiles. Thus the template 
correction functions are generated in a two-step process 
( Appendix lAl summarises these corrections for a particu- 
lar example discussed shortly). The radius R also defines 
the FoV used in determining the overall normalisation ra- 
tio a when estimating event excess and significance (e.g. 
according to Eq. QJ). If the total number of events in 
each w regime out to a radius R is given by P s and Pb 
respectively, the template normalisation a is given by: 

a = Ps/P b (4) 

Alternatively, one may also integrate out to R, the poly- 
nomial functions (un- normalised) p s and pb, to obtain a 
similar result. Typical a values achieved are quite simi- 
lar to that of the displacement background model. The 
template model philosophy is illustrated in Fig. |3] A 
caveat is that the normalisation region must encompass 
all source regions. This will set a natural limit on the 
size of any source region, although systematic effects (sec- 
tion [5J usually dominate for large source sizes. Note that 
the template model correction terms generated from a 
given dataset should strictly be applied only to these same 
data. 

4. Template Model Performance on Real Data 

2D skymaps of excess significance are shown in Fig. 0] 
illustrating the template model performance at various 
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(a) Significance S: Template no correction 
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(c) Significance S: Template full correction 
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Fig. 4. 2D Skymaps and ID distributions of significance S using the template background estimate after various 
correction levels, (a) no correction, Wi — 1.0 (b) radial correction only, /(Az;) = 1.0 (c) full correction as per Eq. 
([SJ. White star: Nominal position (SIMBAD) of Tycho's SNR, White dots: tracking positions, Dashed line: Gaussian 
(/i=0.0, (7=1.0), Solid line: Fitted Gaussian. Events are summed within a circle of radius 0.12° at each grid position 
(120x120 grid), in 0.05°x0.05° steps. 
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Events 




Fig. 3. Illustration of the template background geometry 
with the FoV and normalisation region (defined by radius 
R) in relation to various source positions. The template 
background and '7-ray-like' events are defined by the w 
(MSW) regimes l.3<w < 1.5 and w <1.1 respectively. 
Shown here un-normalised, the template regime contains a 
factor 3~4 more events than that of the 7-ray-like regime. 
Source regions of interest are defined by a cut on some di- 
rectional parameter, e.g. 9 < cut , and are spatially coinci- 
dent in both w regimes. The normalisation region should 
be sufficiently large to encompass all such source regions. 



stages of correction. About 50 hours of CT-System ex- 
posure on Tycho's SNR out to a radius R = 2.24° (VK0°) 
were used in generating these skymaps. Previous analysis 
of these data are described in Aharonian et.al I2001bl At 
each grid point in 0.05° steps, a 8 < 0.12° cut defines both 
the source (which is applicable for point-like sources) and 
template background regions. The significance S is esti- 
mated using Eq. ( ^) . It is apparent the Tycho skymap 
statistics approach that of a Gaussian distribution only 
after correcting for both the radial and zenith-correlated 
systematics. A quite strong zenith-correlated gradient of 
9% per degree is present, which if not corrected, leads 
to biases in the excess significance peaking at the decli- 
nation edges of the FoV (Appendix ^ summarises these 
template correction terms in more detail.). The bin statis- 
tics in these skymaps are somewhat correlated, and lead 
to a slight overestimate in the width a of the ID distribu- 
tion of significances. The level of this overestimate is taken 
from a 2D skymap of significances from bins that have no 
overlap with neighbouring bins. The ID distributions of 
significances for this skymap is a = 0.93 ± 0.03, implying 
that the template background is slightly over-corrected, 
effectively underestimating slightly the significance in the 
FoV. Overall one may conclude however that the template 
model is quite suitable for point source searches over the 
FoV in these data. In this example, the minimum num- 
ber of telescopes required for stereoscopic direction recon- 
struction and w calculation has been chosen at n tc \ > 2, 
and the stereo reconstruction algorithm #1 of Hofmann et 
al. ( 1999) was used. The template model performs equally 
well however under higher choices of n te i (up to 5), vary- 
ing choices of FoV/normalisation radius R, and various 
reconstruction algorithms. 



The template model was also tested with a number 
of a priori chosen positions in the FoV corresponding to 
established 7-ray sources. Presented in Table H are the 
excess significances for a representative sample of TeV 7- 
ray sources with a variety of signal to CR background 
ratios. Results using background estimates from both the 
template and displacement background models are given 
for comparison. In all cases a normalisation radius for the 
template background model of R — 2.24° and minimum 
number of images ritd > 2 have been used. The tem- 
plate model provides derived excesses s — ab quite con- 
sistent with those obtained from the displacement model. 
Since the example sources of Table ^ lie quite close to 
their respective tracking positions, the zenith-correlated 
correction term of the template background has little ef- 
fect, and hence results only for the fully corrected tem- 
plate background are shown. For the Crab, a slight mod- 
ification regarding the template correction is made to the 
fit limits of the polynomials p 8 and pt, to exclude re- 
gions containing this strong source. The polynomial fits 
are affected by strong sources near (< 0.5°) the track- 
ing position and these should be avoided. It should be 
stressed that the derived excesses and significances for 
the Cas A and H1426+428 and 1ES1959+650 here are 
generally lower then those of presently published analyses 
(Aharonian et al. I2001al l2T)Tl2l Horns et al. EES), which 
employed the more sensitive reconstruction algorithms #2 
and #3 (Hofmann et al. 1999), and in some cases a higher 
minimum n to i- The template model does reproduce pub- 
lished results when employing the corresponding analyses. 

Finally, a particularly useful feature of the template 
model is demonstrated against a problematic systematic 
effect arising when camera pixels view bright stars in the 
FoV. Increased noise fluctuations (and corresponding an- 
ode currents) in camera pixels viewing bright stars can 
distort derived image parameters. Affected pixels, usu- 
ally the same for all five cameras, are therefore removed 
from a camera trigger dynamically during data taking 
and also from image analysis at a software level. Any 
pixel with anode current > 3fiA is not considered in es- 
timation of image parameters and consequent event di- 
rection reconstruction (Bulian et al. 119981 Daum et al. 
1997). A residual effect of this latter step is to leave a 
systematic deficit of reconstructed event directions in a 
region broadly centred on the offending star. This effect 
was verified by artificially removing certain pixels before 
image paramctrisation. In practise, up to two to three 
pixels may be removed during data taking at any given 
time depending on the star brightness. The Crab field, 
with C-Taurii (mB=2.9) and the Cas-A field with AR- 
Cas (me =4.8), are ideal examples to illustrate such star 
deficits. Quantitatively, star deficits amount to a reduc- 
tion in CR events up to ^20% over a diameter ~ 0.2°. 
Table summarises the template (after full correction 
only) and displacement background model performance 
at the £-Taurii and AR-Cas positions. Since the template 
background closely follows the CR background response 
of gamma-ray-like events, star deficits are considerably re- 
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Source 


RA [hr] Q 


Dec [deg] a 


cut [deg] 


s b 


Q 


s — ab 


S [*]" 




— Template background (full correction: Eq ( 








Cas-A 


no onri 
zo.oyu 


58.82 


0.12 


1364 4649 


0.258 


165 


+4.2 


H1426+428 


1 A A7f\ 


42.67 


0.12 


385 954 


0.295 


104 


+5.1 


1ES1959+650 


20.000 


65.15 


0.12 


571 2236 


0.210 


101 


+4.1 


Crab 


5.576 


22.01 


0.12 


1756 884 


0.282 


1507 


+48.9 






— Displacement background — 








Cas-A 






0.12 


1364 6092 


0.197 


164 


+4.2 


H1426+428 






0.12 


385 997 


0.283 


103 


+5.1 


1ES1959+650 






0.12 


571 2437 


0.201 


81 


+3.2 


Crab 






0.12 


1756 1281 


0.226 


1466 


+48.2 



a. J2000 epoch. 

b. Statistical significance from Eq. f[TV 



Table 1. Event excesses s — ab and significances S derived using the template and displacement background models 
for a number of established TeV 7-ray sources. The source region is defined by the cut 6 < C ut- s: 7-ray-like events; 
b: background estimate; a: normalisation. Datasets are from the HEGRA CT-System archive; Cas A (212 hours), 
H1426+428 (42 hours, 1999/2000), 1ES1959+650 (94 hours, 2000/2001), Crab (33 hours). 



Star 



RA [hrs] a Dec [deg] a 



ut [deg] 



ab S \a} c 



Template background (full correction: Eq. (EJ) 

C-Taurii 5.627 21.14 2.4 0.15 970 3399 0.290 -16 -0.4 

AR-Cas 23.501 58.55 4.8 0.15 1294 5293 0.258 -72 -1.7 



C-Taurii 
AR-Cas 



Displacement background 

2.4 0.15 970 6241 0.201 -290 -7.6 
4.8 0.15 1294 7358 0.213 -273 -6.4 



a. J2000 epoch. 

b. B-band star magnitude. 

c. Statistical significance from Eq. f[Tl. 



Table 2. Event excesses s — ab and significances S at star positions within the Crab (£-Taurii, 161 hrs) and Cas-A 
(AR-Cas, 212 hrs) fields of view. Compared are results using the template and displacement background models. 



duced when using the template model. The displacement 
background however, which does not include the star re- 
gion, produces a strong negative bias. Conversely, if dis- 
placement background regions (applied to another source 
region) were to overlap star deficit (s), some positive biases 
could be introduced. A large number of bright stars in 
the FoV would in general render application of displace- 
ment background models problematic. Overall the tem- 
plate model appears able to match localised systematic 
changes in response of 7-ray-like CR events in the FoV. 

5. Limitations 

The template model as designed here employs corrections 
only for large-scale systematic changes in CR response, 
correlating radially from the tracking position and with 
event zenith angle. Higher-order systematics (correlating 
with alternate system parameters) arising from, say, differ- 
ences in telescope camera response can also reasonably be 
expected, despite long-term efforts to maintain a homoge- 



neous response in the HEGRA CT-System (see Puhlhofcr 
et al. I2UU3|I . Since systematic effects will grow in propor- 
tion to the size of the source region and also exposure, 
these are better-revealed in deep exposure datasets. For 
the Tycho's SNR dataset however, significance skymaps 
generated with larger source sizes in mind (e.g. 6 cut up to 
0.6°) do not show behaviour significantly differing from 
a Gaussian distribution. An upper limit to the system- 
atic error in the derived excess s — ab across the Tycho 
FoV can be estimated by using subjectively-chosen 'hot' 
and 'cold' regions (of size ~ 1.0°) of the fully corrected 
skymap. In this way a systematic error of ~ ±4% in units 
of the normalised background ab is calculated. 

The Tycho's SNR dataset is representative of many in 
the CT-System archive of order 50 hrs exposure, with >4 
telescopes of the CT-System in operation. Such an expo- 
sure is generally found not sufficient to reveal higher-order 
systematic effects. The Cas A dataset comprises a deeper 
exposure (> 200 hrs). The fully corrected template model 
skymap of these data (Fig EJ shows a strong deviation 
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Fig. 5. Skymap (2D & ID) of significance S using the fully corrected template background estimate for Cas-A data. 
White star: Nominal position (SIMBAD) of Cas-A, White dots: tracking positions. Dashed line: Gaussian (/j,=0.0, 
cr=1.0). Details for skymap construction are the same as for Fig0J 



from Gaussian behaviour (ID a = 1.126 ± 0.008) already 
when searching for point sources (0 cut < 0.12°). A 'hotter' 
group of bins is seen at RA> 23.5 h, Dec 58° to 59°, rep- 
resenting an additional systematic excess of order few per- 
cent in units of normalised background. The source of this 
systematic is as yet not fully understood. It may be linked 
to a 'missing' telescope (CT2 is not present in large frac- 
tion of these data) which is uncovering inhomogeneities 
in trigger response correlating with certain telescopes, al- 
though this gradient does not align with the CT-System 
geometry. It is seen in these data that artificial removal 
of other telescopes for example, can also introduce fur- 
ther similar-scale effects which may or may not be related. 
Changes in average skynoise across the FoV may also con- 
tribute. Care should therefore be taken when searching for 
sources in affected regions, employing alternative back- 
ground estimates such as displacement backgrounds as a 
check. All background methods will of course also be af- 
fected by large-scale systematics of various types and so 
the above-mentioned problem is not particular to the tem- 
plate model. Extra correction terms could be added to the 
template model to accommodate such higher order effects 
but this has not yet been explored in detail. Experience 
has shown that for large exposures, higher order system- 
atic effects can be averaged out or considerably reduced 
by employing a variety of tracking positions, for exam- 
ple in a mosaic fashion. All of the example datasets pre- 
sented here employed the so-called wobble mode in which 
only two tracking positions were alternated ±0.5° in dec- 
lination with respect to the primary source of interest. 
Investigation of higher order effects and their origin is 
presently continuing, in the context of a search for TeV 



sources from the entire CT-System archive (Puhlhofer et 
aLimH Aharonian et al. l2TM|> . 

The estimation of energy spectra from source excesses 
requires that the background estimate be comprised of 
events with a similar distribution of estimated energies 
to that of events from the source region. The template 
background estimate however does not meet this criterion 
since the w values comprising the background estimate 
will differ in estimated energy from those in the 7-ray-like 
w regime. Thus alternative background models should be 
used for spectral analysis. 

6. Conclusions 

Described in this paper is a new method to estimate the 
CR background in 2D searches for sources of TeV 7- 
ray emission with ground-based detectors. This so-called 
template background estimate, demonstrated here with 
HEGRA CT-System data, employs a subset of CR events 
normally rejected according to the Cerenkov image shape 
parameter mean-scaled-width w . These template events 
spatially and temporally overlap with CR events consid- 
ered 'gamma-ray-like', and thus no dedicated OFF source 
observations are required. Applying the template model 
successfully to HEGRA CT-System data involves correc- 
tion over of the field of view for: (1) differences in radial 
response and (2) differences in a zenith-correlated gradi- 
ent, both between the CR events of the two w regimes. 
A particularly useful feature of the template model is its 
ability to compensate well for localised systematic changes 
in CR event density due to the presence of stars in the field 
of view. It is well suited to cases where many stars and 
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TeV sources are present. The template model applicabil- 
ity is limited by the presence of systematic gradients in 
CR event density which are presently not corrected. All 
background models will however suffer to various extent 
from weak, large scale systematics. Further investigation 
of these aspects is ongoing. Systematic uncertainties in 
the derived event excess of less than 4% of the normalised 
background are achieved when searching for sources of size 
less than ~ 0.6° radius in many datasets. 

It should be possible to extend the template model 
philosophy to analyses which use other means to reject 
CR background, for example maximum likelihood meth- 
ods, multi-dimensional cluster analyses, applied to either 
single or multi-telescope systems. Issues relating to sys- 
tematic gradients over the FoV and star deficits will be 
important for the next generation IACTs such as H.E.S.S., 
VERITAS, CANGAROO III and MAGIC. In these new 
systems, CR event rates of order 500 Hz are expected, 
likely revealing higher-order systematic gradients in re- 
sponse after observation times significantly shorter than 
that required in HEGRA CT-System data. Furthermore, 
the larger mirror areas of these future experiments will 
also increase considerably the number of stars that could 
cause detrimental systematics in the FoV response. The 
template background model could therefore be quite use- 
ful in 2D skymap construction, in future searches for TeV 
7-ray sources. 
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Appendix A: Template Background Corrections 

Summarised here are the two template background cor- 
rection terms calculated specifically for the HEGRA CT- 
System Tycho's SNR dataset (Fig.©. 

1. Radial Correction: The radial event density for the 
entire dataset, centred on the various tracking posi- 
tions may be fitted with 8th-order polynomials out to 
a radius R — 2.24°: (1) p s for the 'gamma-ray-like' 
(w <1.1), and (2) pt for the template regime (1.3<u> 
<1.5): 

8 

P.,& = X> (A.l) 

i=0 

where for Tycho's SNR data one fits: 



p s pt 



Po 


215.3 ±0.7 


815.3 ±0.9 


pi 


17.0 ±0.8 


160.8 ±1.9 


P2 


-148.7 ±0.1 


-606.2 ± 0.8 


P'i 


148.00 ± 0.02 


576.80 ±0.19 


Pi 


-62.340±O(10 -3 ) 


-243.60 ± 0.04 


P5 


7.987±O(10" 3 ) 


31.95 ±0.01 


P6 


1.661±O(10 -4 ) 


6.117±O(10 -3 ) 


P7 


-0.5548±O(10 -5 ) 


-2.103±O(10~ 4 


PS 


0.0418±O(10~ 4 ) 


0.158±O(10~ 3 ) 



2. Zenith-Correlated Correction: A linear function 
is used to describe the ratio of event densities, as a 
function of the zenith angle difference (Az) between 
the event and tracking positions, of the gamma-ray- 
like and template w regimes,: 

/(Az) = (0.998 ±0.006) ± (0.898 ± 0.006) Az (A.2) 

A linear gradient of order 9% per degree along the 
zenith axis is therefore present in the Tycho's SNR 
data. Separately, one can also fit functions to each id 
regime; f(Az) s and /(Az)t,, and note that the function 
/(Az) = /(Az) s //(Az) b : 

f(Az) s = (0.999 ± 0.006) + (0.107 ± 0.006) Az 
f(Az) b = (1.000 ± 0.006) ± (0.015 ±0.006) Az 

The zenith-correlated gradient therefore arises pre- 
dominantly from CR events in the gamma-ray-like w 
regime. 

This gradient actually comprises two components: (1) 
the naturally expected event rate dependence upon 
zenith angle for a Cerenkov detection instrument, and 
(2) using in the construction of the w parameter 
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which employs a zenith-dependent scaling, the track- 
ing zenith angle instead of that for each individual 
event. Here, the template model is demonstrated for 
a type of worst-case-scenario, without removing com- 
ponent (2). Component (1) is usually of order few per- 
cent per degree for data taken at zenith angles such in 
the Tycho's SNR dataset. 

It should also be noted that Eq. ( ^) takes no account 
of statistical errors in the normalisation a, and the tem- 
plate correction terms. These errors are however generally 
less than one percent (for obs. time > 10 h) and can be 
neglected given the high significance levels to which these 
parameters are estimated. 



